The Effects of Molecular Anions on the Chemistry of Dark Clouds 

Catherine Walsh 1 , Nanase Harada 2 , Eric Herbst 3 and T. J. Millar 1 

cwalshl3@qub .ac.uk 

ABSTRACT 

We have investigated the role of molecular anion chemistry in pseudo-time-dependent 
chemical models of dark clouds. With oxygen-rich elemental abundances, the addition 
of anions results in a slight improvement in the overall agreement between model results 
and observations of molecular abundances in TMC-1 (CP). More importantly, with the 
inclusion of anions, we see an enhanced production efficiency of unsaturated carbon- 
chain neutral molecules, especially in the longer members of the families C n H, C n H2, 
and HC n N. The use of carbon-rich elemental abundances in models of TMC-1 (CP) 
with anion chemistry worsens the agreement with observations obtained in the absence 
of anions. 

Subject headings: astrochemistry — ISM:abundances — ISM:clouds — ISM:molecules 



1. INTRODUCTION 

The possible existence of molecular anion s ( negatively ch arged ions) in the in ter stellar medium 



was discussed by iDakarno fc McCravi fll973]), ISarrd fll980|), and iHerbstJ (Il98lh . iHerbstl (I198lh 
calculated that the large electron affinity of carbon-chain molecules and hydrocarbon radicals would 
lead to efficient radiative electron attachment rate coefficients for species with more than 4-5 atoms 
resulting in anion-to-neutral ratios on the order of a few percent. Confirmation of this hypothesis 
was achieved followin g the measurement of th e rotational spectrum of the molecular anion, CgH~ , 
in the laboratory by iMcCarthv et al.l (120061) . This measurement allowed the verification, from 
existing astronomical spectra (jKawaguchi et al.l 1 19951 ) , of the presence of CqH~ in the envelope of 
the carbon rich evolved star, IRC+10216, at an abundance l%-5% that of the neut ral molecule. A 



successful astronomical search for CqH in TMC-l(CP) bv lMcCarthv et al.l (|2006l ) determined an 
anion-to-neutral ratio of ~ 2.5% in this source. 
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The rotational spectra of C4H , CsH , and C3N have since been measured in the laboratory 



(Gu pta et al 



2007 



Thaddc us et al.ll2008l ) wit h subsequent detectio ns of C4H (jCernicharo et al 



2003), C 8 H™ (IRemiian et all2007h. and C^N™ (jThaddeus et all200a ) in the envelope of IRC+10216 



and CsH in TMC-1 (jBriinken et al.l 120071 ), Prompted by the observation of various c arbon 



chain molecules with large abundances in L1527, a class O/I protos tar. by ISakai et al.1 (|2008l ) 



CgH™ was searched for and confir med to be present i n this source (jSakai et al.l 1200^ ) followed 
clo sely by the d e tectio n of C4H™ (jAgundez et al.1 120081 ) . A survey of galactic molecular sources 
by iGupta et al.1 (|2009l ) detected CeH™ in two futher sources, the pre-stellar cloud, L1544, and 



the protostell ar object, L1521F , sugg estive of the ubiquitousness of molecular anions and CqH~ 
in particular. ICernicharo et al.l (|2008l ) have attributed a series of rotational lines observed in the 
envelope of IRC+10216 to a new molecular anion, C5N - , although confirmation of this detection 
awaits successful measurement of the rotational spectrum of the molecule in the laboratory. 

Foll owing the id e ntifica tion of CgH™ in the envelope of IRC+10216 and the detection in 



TMC-1, Millar et al 



(2007) expanded the current release of the UMIST Database for Astro- 



chemistry (IWoodall et al.l 120071 ) and the existing chemical network for the envelope of IRC+10216 



dMillar et alJbood ) ;o include the molecular anions, C n with n = 5 - 10 and C n H with n = 4 



10, 

using calculated rate coefficients for radiative electron attachment. Their results were limited to 
calculated anionic and precursor neutral column densities. For TMC-1 (CP), they found that the 
cal culated 'early time' Cr H™ column density of 1.35 x 10 11 cm -2 was very close to that measured 



bv lMcCarthv et al.l (|2006l ). 1 x 10 cm™ . The calculated anion-t q-neutral ratio of 5 . 2% wa s also in 
good agreement with the observed value of 2.5%. More recently, lHarada &: Herbsti (|2008l ) investi- 
gated the chemistry of the protostellar ob j ect, L 1527, concentrating on the formation and depletion 
of molecular anions. IWakelam &: Herbsti (|2008l ) considered the effects of negatively charged PAHs 
(polycyclic aromatic hydrocarbons) on dense cloud chemistry and found that PAH™ can replace 
electrons as the dominant negative charge carrier, an effect which depends strongly on the size of 
the PAHs considered. 

The presence of molecular anions at significant abundances in the sources discussed above 
suggests that anions play a more significant role in the chemistry of astrophysical environments 
than previously believed. In this paper, we investigate the effects that the inclusion of molecular 
anion chemistry has on the chemistry of dark clouds by running simple pseudo-t ime-dependent 



mode ls using the most recent release of the UMIST Database for Astrochemistry (jWoodall et al 



20071 ). henceforth referred to as 'Rate06', and the latest OSU (Ohio State University) network, 



osu.03.2008 (see http ://www.physics.ohio-state.edu/~eric/[ ). We compare our calculated fractional 
abundances for models both excluding and including anions with fractional abundances observed 
towards the cyanopolyyne (CP) peak of the dark cloud core, TMC-1, to determine the role that 
these species play in the gas-phase chemistry of dark clouds. 

In Section [21 we describe our anionic chemistry and the initial chemical abundances and phys- 
ical conditions adopted for our standard model. Section [3] contains detailed results obtained with 
and without anions for both networks, discusses the effects of anions on the chemistry of selected 
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classes of neutral species and compares our results with observations towards the molecular-rich re- 
gion, TMC-1 (CP). Section U] shows how our results are affected by the use of carbon-rich elemental 
abundances. Finally, in Section [5] we summarize our main findings. 



2. CHEMICAL MODEL 



We have expanded both Rate06 and the OSU network to include the molecular a nions C n (n 
3 - 10), C n H~ (n = 4 - 10), and C n N~ (n = 3, 5), building upon the previous work bv lMillar et al 



(|2007l ) by increasing the number of molecular anions studied and expanding the reaction set for 
each anion to inc orporate new reaction ch annels and rate coefficients from recently published ex- 
perimental work (jEichelberger et al.l 120071 ) . We use the dipole-enhanced version of Rate06, which 
includes an enhancement of ion-neutral rate coefficients at low temperatures for reactions where the 
neutral has a large, permanent dipole mo ment (> 0.9 Debye). The resulting rate coeffi cients have 



a T 1 / 2 dependence at low temperatures ( Herbst h Leunall986l ; IWoodall et al.ll2007l ). A similar 



but not i dentical enhancement for the rate coefficients of ion-dipolar reactions is used in the OSU 
network (jHerbst LeungHl98a ). 



The anions are formed primarily via radiative electron attachment, 



X + e" 



X + hv, 



(1) 



with u pdated attachment rate coefficients from recent theoretical calculations bv lHerbst fc Osamura 
(|2008l ). The hydrocarbon anions are destroyed through reactions with H, C, O, N, by photo- 
detachment, and by mutual neutralization with abundant cations such as C + . The destruction 
channels are listed in reactions ([2]) to ([7]) with m = or 1: 



CnHj, 

C n H E 
C n H n 
C n H n 



+ H 

+ C 
+ 
+ N 



C n H m 



" + hv 
+ X+ 



C n H m+ i + e , 
C n +iH m + e~, 
CO + C n _i H m , 
CN + C n _iH m -, 
CN" + C n _iH m , 
C 3 N- + C n _ 3 H m , 
C5N + C n _5H m , 
H m C n N + e-, 
C n H m + e , 
C n H m + X. 



(2) 
(3) 
(4) 
(5) 



(6) 
(7) 



We have taken the ra te coefficients for reaction ([2]) from iBarkholtz et al.l (|200ll ) and for reactions 
([3]), d3|), and ([5]) from lEichelberger et al.l (|2007l ) for the species for which they were measured. In 
the absence of laboratory data, we have used estimated rate coefficients based on the experimental 
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measuremen ts. The rate coeffic i ents a nd reaction channels for reaction ([5]) are particularly uncer- 
tain because lEichelberger et al.l (|2007l ) were unable to detect associative detachment as a reaction 
channel, a potentially important synthetic route to the HC n N molecules. We discuss the implica- 
tions of these uncertainties in ^3.2.51 The photo -detachment rate coefficients (reaction [6] ) adopted 
are those listed in Table 1 of lMillar et al.l (|2007l ). although in our model, this is an unimportant de- 
struction mechanism as TMC-1 is well-shielded from external sources of UV radiation. The mutual 
neutralization reac tions (reaction ( 7 ) ) each have a rate coefficient k = 7.5 x 10 -8 (T/300) -0 ' 5 cm 3 
s _1 as included by lHarada Sz Herbstl ( 20081 ) in their chemical model of L1527. We have assumed 
in all instances where appropriate that the products are obtained by a simple electron tra nsfer. 
The CaN" chemistry we have added includes that listed in Table 1 of lHarada fc Herbstl (|2008l ) plus 
formation routes via the molecular anions , C n ~ and C n H~ (eq. [5]), as observed in the experi- 
mental results of lEichelberger et al.l (|2007l ). We have included C5N - using the same destruction 
reactions and rate coefficients as for CsN - with a formation reaction r ate coefficient (react ion [TJ) 



calculated by E. Herbst (private communication) using the method of lHerbst &: Osamura 
to be 1.25 x 10- 7 (T/300)-°- 5 cm s . Due to the inclusion of C10 and C10H , we extrapolated 
both networks to include further hydrocarbon species containing 10 carbon atoms hence restricting 
the hydrocarbon chemistry to molecules with 11 carbon atoms or less. 

We model the gas-phase chemistry of a dark cloud core by treating the source as an homoge- 
neous, isotropic cloud with the constant physical parameters, n(H2) = 10 4 cm -3 , T = 10 K, and 
A v = 10 mag, using a cosmic-ray ionisation rate of 1.3 x 10~ 17 s~ 1 , and performing pseudo-time- 
dependent calculations, in which the chemistry evolves from initial abundances that are atomic 
except for molecular hydrogen, which is produced efficiently at an earlier stage. We employ the 
commonl y used oxygen- r ich low -metal elemental abundances, as listed in Tableland suggested orig- 
inally by iGraedel et al.l (119821 ). Because pseudo-time-dependent gas-phase models of dark clouds 
are relatively predictive compared with more complex physical models, using such a simplified 
model helps us to accurately identify particular reactions or, more commonly, systems of reactions 
that are important and require further study. 

Any comparison between our modelled results and observations in TMC-1 (CP) could, of 
course, be compromised by our use of a simplified model for a dark cloud. We would expect 
that employing a more detailed physical model could lead to different results from those presented 
here; however, such mod els are certainly not unique and th erefore bring their own limitations. For 
example, recent work bv lHassel et al.l (|2009. in preparation! ), in which a chemical dynamical- model 
following both the gas-phase and surface chemistry as a dense cloud forms via a shock, finds little 
difference in the gas-phase abundances of large molecules although the time-scales for formation 
are certainly different. Chemical factors, such as our neglect of accretion of gas-phase species onto 
dust grains, will also affect gas-phase abundances, although the timescales we are discussing here, 
~ 10 5 yr, are only just at the stage where accretion becomes important for the conditions we have 
adopted. 
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3. RESULTS 

Table [2] shows our calculated abundances using both networks for an assortment of molecules 
at the times of best agreement (typically 1 — 2 x 10 5 yr) for TMC-1 (CP) along with the observed 
values or upper limits. Model results are shown for calculations with and without anions. As can be 
seen, there is little difference between the results obtained with Rate06 and with the OSU network. 
We first discuss the anions and then consider families of neutral species strongly affected by the 
inclusion of anions. Throughout this section, we display the results from the models using Rate06 
only because the results from the models using the OSU network are similar. 



3.1. Anions 

The calculated fractional abundances of the molecular anions, C n H~ for n = 4, 6, 8, and 10, 
C n ~ for n = 4, 6, 8, and 10, and C n N _ for n = 3 and 5, along with the neutral precursor or each 
anion, are plotted as a function of time in Figures [T] (a), [1] (b) and [I] (c), respectively. Because the 
main route of formation of the anions at all times is radiative electron attachment, the abundances 
of the anions follow that of their neutral precursors with both generally inversely proportional 
to molecular size. An exception is the abundance of C4H~~, which at ~ 10 5 yr is comparable 
to that of Cf;H~ becau se of the slow radiative electron attachment rate coefficient calculated for 



C4H (IHerbst &; Osamurall2008l ). The larger abundance of C5N - compared with C3N - at all times, 
despite the smaller abundance of C5N compared with C3N, is due to the much larger electron 
attachment rate coefficient calculated for C5N (~ 10~ 7 cm 3 s _1 compared with ~ 10~ 10 cm 3 s _1 ). 
Although the reaction rate coefficient for the formation of CsN - via dissociative attachment of 
HNC 3 is two orders of magn itude larger than that of the radiative electron attachment (see Table 



1 of lHarada & Herbstl (|200J)) the abundance of HNC 3 only reaches ~ 0.01 that of C 3 N after ~ 10 5 
yr when dissociative attachment becomes comparable with radiative attachment. The abundances 
of both the anions and neutrals peak around 10 5 yr and the destruction of the anions is dominated 
by reactions with H, C, O and C + (reactions [2j, [3], [4], and [Zj, respectively). 

In the absence of anions, the main destruction channels for C n H are via reactions with C, O, 
and N: 

C + C n H -» C n+ i+H, (8) 
+ C n H -> C n _iH + CO, (9) 
N + C n H -> C n _iH + CN. (10) 

The inclusion of anions introduces radiative electron attachment (reaction pQ) as a destruction 
mechanism and it is this channel that dominates the destruction of the C n H molecules at all times, 
excepting C4H for which the destruction via the reaction with O is comparable. 



Anionic chemistry also introduces new formation routes for the C n H radicals. Ordinarily, these 
radicals are formed primarily via dissociative electron attachment to hydrocarbon cations as well 
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as via carbon insertion reactions (in reaction pT] m = 2 or 3): 

C n H m + +e- -> CnH + H^!, (11) 
C n+ iH + + e- -» C n H + C, (12) 
C + C n _iH 2 -» C n H + H. (13) 

New formation routes become important past ~ 10 3 yr when appreciable abundances of molecular 
anions have built up and C n H is formed via associative electron detachment from carbon-chain 
anion and hydrocarbon anion precursors (see reactions [2] and [3]): 

H + C n - -» C n H + e", (14) 
C + C n _iH- -» C n H + e-. (15) 



Radiative electron attachment of the bare carbon-chain neutrals, C n , as a destruction mecha- 
nism, only dominates prior to 10 4 yr, at which point reactions with O become comparable to this 
channel. As with the hydrocarbon radicals, the inclusion of bare carbon-chain anions opens up an 
additional formation route to the bare carbon-chain neutrals via associative electron detachment 
(reaction [3]): 



c„_r + c 



C n + e" 



(16) 



As can be seen in Table [2j the abundances for CqH~ and CsH - are over-predicted at the time 
of best overall agreement by around an order of magnitude, with the abu ndance for C4H" around 
50 times greater than the estimated upper limit by lThaddeus et al.1 (120081 ). There is a sharp decline 
in the predicted C n H~ abundances at somewhat later times, however, as can be seen in Figure [J (a). 
The abundance of C4H - only falls below the upper limit estimated by lThaddeus et al.1 (|2008l ) after 
~ 10 6 yr indicating that our electron attachment rate for C4H is too large or that we are neglecting 
further destruction mechanisms for this species. In contrast to the anions already discussed, the 
predicted abundanc e of CaN~ at the t ime o f best overall agreement is in accordance with the upper 
limit determined by iThaddeus et al.l (|2008l ) . 



The calculated anion-to-neutral ratio, C n H~/C n H for n = 4, 6, 8, and 10 is plotted as a 
function of time in Figure along with the observed ratios for CqH and CsH. The calculated 
ratio for all species reaches a minimum around 10 5 yr, the same time that the peak in individual 
fractional abundances is reached. The model results for Rate06 yield 5.2% for the CgH ratio and 
3.9% for the CsH ratio. These compare with the observed values for CgH and CsH of 1.6% and 
4.6%, respectively. 



3.2. Effects on Other Species 



The main aim of this paper is to investigate the effects that the inclusion of anions has on the 
abundances and evolutionary behavior of other species in established chemical models. We have 
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identified several families of species that are affected; these include the neutral precursors of the 
anions, C n , C n H, and C n N as well as the families of molecules C n H2 and HC n N. The sensitivity of 
species in these families to the inclusion of anions can be seen in Table [2] for both networks. 

3.2.1. C n 

Figure [3] (a) displays the fractional abundances of the neutral C n molecules for n = 4, 6, 8, and 
10 as a function of time for the models with and without anions. It is immediately apparent that 
the abundances of Cq, Cg, and Cio are enhanced at all times. The relative enhancement is more 
significant with increasing molecular size with little change in the abundance of C4, in contrast to 
that of Cio, which increases by around 2 orders of magnitude in the presence of anions (see Table 
|2J). As discussed previously in §3 . 1 1 introducing anions to the chemistry includes a new formation 
route via associative electron detachment (reaction [16] ) whereby larger neutral carbon-chains can 
build up from smaller anionic precursors through repeated insertion of a carbon atom. 

3.2.2. C n H 

Figure [3] (b) shows the fractional abundances of the C n H molecules for the models with and 
without anions. As in Figure[3](a) for C n , there is an obvious enhancement of the larger hydrocarbon 
radicals, CgH, CsH, and C10H, at all times prior to 10 6 yr with the enhancement increasing with 
molecular size. The addition of anions introduces two new formation mechanisms via associative 
detachment from bare carbon-chain anions (reaction p3]) and hydrocarbon anions (reaction |15|). 
again introducing a new route to the formation of larger hydrocarbon radicals from smaller ones. 
At the time of overall best agreement, as seen in Table El the inclusion of anions enhances the 
abundance of C@H from a value which is ~ 3 times smaller than observed to one ~ 3 times larger 
than observed, while the abundance of CsH is enhanced from a value ~ 4 times smaller to one as 
14 times larger than observed. There is no enhancement of the C4H abundance, which remains a 
factor of ~ 5 below the observed value. As in the case for the hydrocarbon anions, the abundances 
of both CqR and CsH reach their best agreement with observation at slightly later times, while 
that of C4H diminishes further. 

3.2.3. C„N 

The fractional abundances of the C n N radicals, with n = 3, 5, 7, and 9, are plotted in Fig- 
ure [3] (c), as a function of time. The abundances of all of these species are enhanced in the model 
with anions, compared with the results for the model without anions, prior to ~ 3 x 10 5 yr. As 
with the carbon chains and hydrocarbons, the biggest relative enhancement is seen in the larger 
species, C7N and C9N, with the peak fractional abundance of C9N increasing by around 2 orders of 
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magnitude from ~ 10 to ~ 10 . When anions are not included in the model, the main routes 
to the formation of C n N for n = 5, 7, and 9 are 

H 2 C n N+ + e- -» C n N + H 2 , (17) 
H 3 C n N+ + e- -► C n N + H 2 + H, (18) 
N + C n H -> C n N + H, (19) 

with destruction via reactions with O and N atoms. Reaction (|17|) dominates the formation up 
to ~ 10 5 yr after which reaction (|19p is comparable to or dominates over reaction (|17p . For C3N, 
reaction (|19p is the main route of formation at all times. 

The inclusion of anions introduces a new formation route via the bare carbon chain anions, 

c„- 

N + C n " -» C n N + e-, (20) 

and it is this route that dominates the formation for n = 5, 7, and 9 up to ~ 10 5 yr after which 
reactions ()17p and (|19p take over. As in the case without anions, the main formation route for C3N 
is through reaction (|19p . As the abundances of the C n H radicals are enhanced with the inclusion 
of anions this also increases the influence of reaction (|19p over that of reaction (|17|) in contrast to 
the case without anions. 



3.2.4- C« H2 

Whilst the effect on the abundances of C n and C n H due to the inclusion of their respective 
anions may come as little surprise, other families of species are also affected. Figure 0] (a) shows 
the time-dependent fractional abundance of C n H 2 for n = 4, 6, 8, and 10. The abundances of all 
species are enhanced at all times in the model with anions included, with a more marked effect for 
CsH 2 and CioH 2 . In the absence of anions, the C n H 2 species form via dissociative recombination, 
e.g. reactions (|2ip and (|22p . or via neutral-neutral reactions, 

C n H 3 + + e- - C n H 2 + H, (21) 
C n H 4 + + e- -» C n H 2 + H 2 , (22) 
C n _ 2 H 2 + C 2 H -» C n H 2 + H. (23) 

Adding anions to the chemistry opens a new synthetic channel via associative detachment of the 
hydrocarbon anions, C n H~, with H atoms (reaction [2]) and it is this channel which dominates 
up to ~ 10 6 yr. The C n H 2 species in the dominant polyyne isomer, HC n H, are not observable 
via rotational spectroscopy because this isomer has zero dipole moment. The cumulene carbene 
forms H 2 C n are detectable but in general cannot be formed by associative detachment because the 



reactions are endothermic (IHerbst Osamurall2008l ). 
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3.2.5. HC n N 

Figured] (b) shows the abundances of the HC n N molecules for n = 3, 5, 7, and 9 and again we see 
a marked enhancement in the abundances of all species when anions are included in the networks. 
The synthesis of the larger cyanopolyyne molecules proceeds (as for many neutral molecules) via 
the dissociative recombination of a cationic molecular precursor: 

H 2 C n N+ + e- -» HC n N + H, (24) 
H 3 C n N + + e - -» HC n N + 2H, (25) 

with a secondary route through the reaction 

C n _iH 2 + CN -» HC n N + H. (26) 

The addition of anions to the chemistry enhances the abundances of the C n H2 species which in- 
creases the effect of reaction (I26p . There is also a more subtle influence on the abundances of 
H 2 C n N + and H3C n N + which are produced via a series of different reaction pathways originat- 
ing from hydrocarbon cations, C n H m + , the abundances of which are enhanced in the presence of 
molecular anions. Examples of this type of reaction are 

C 4 H 2 + +HC n _ 4 N -» H 3 C n N + + ^, (27) 
C 2 H 2 ++HC n _ 2 N -» H 3 C„N+ + /u/. (28) 

In addition, the inclusion of C3N - and C5N - introduces new formation mechanisms for HC3N and 
HC5N via associative electron detachment of each anion with H atoms. This reaction has a fast 
rate coefficient, on the order of 10 -9 cm 3 s _1 , which has a large effect on the abundance of HC5N, 
in particular. 

The destruction of the hydrocarbon anions C n H~ via reactions with atomic N (reaction [5]) 
has an associative detachment channel, which leads to the cyanopolyyne molecules for n = 3, 5, 
7, and 9. The rat e coeffi cients for these reactions have been derived from experimental data by 



Eichelberger et al.l (|2007l ) who found that the reactions of anions with an even number of carbon 
atoms proceed slowly, with rate coefficients on the order of 1CP 11 cm 3 s . Branching ratios for 
the products were not determined, although contact with the authors led to their agreement that 
the associative detachment channel to form the cyanopolyynes could have a branching fraction 
of ~ 0.5. For reactions with odd numbers of carbons in the hydrocarbon anions, i.e. n = 3, 5, 
7, and 9, we have assumed that the reactions proceed with a rate coefficient the same order of 
magnitude, with an associative detachment channel that has a branching fraction 50% of the total 
rate. Inclusion of this reaction channel has little effect on the HC n N abundances due to the small 
rate coefficient used; however, we have also run a model with this reaction channel proceeding using 
a larger rate coefficient of 10 -10 cm 3 s _1 , on the basis that the reactions between C n H~ species 
with odd values of n and atomic nitrogen are spin allowed. The results from this "third" model are 
also shown in Figure [4] (b) (dotted lines) in addition to the results from the model without anions 
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(solid lines) and our standard model with anions (dotted-dashed lines). The associative detachment 
channel futher enhances the abundances of HC7N and HC9N in particular, two molecules for which 
existing chemical networks have difficulty reproducing the observed abundances. At the time of 
best agreement in the third model (1.6 x 10 5 yr), the abundances are in vastly better agreement 
with the observed abundances listed in Table [2] although they still lie somewhat below observational 
values. 



3.3. Comparison of Model Results with Observed Abundances in TMC-1 (CP) 



In Table [21 the abundances of species that agree to within one order of magnitude at the time of 
best overall agreement are in normal type, those predicted to be more than one order of magnitude 
too large are in bold italic type and those predicted to be more than one order of magnitude too 
small are in Roman bold type. It is immediately apparent that the inclusion of anions in both 
networks does not have a dramatic effect on the agreement with observation. Specifically, removing 
the anions from Rate06 reduces the percentage of molecular abundances that lie within one order 
of magnitude of those observed from ~ 66% (35 out of 53 species) to ~ 63% (32 out of 51 species). 
The improvement when anions are included is minor but noticeable. 

Important molecules for which the inclusion of anions improves agreement with observation 
are HC5N, HC7N and HC9N (as discussed in detail in §3.2.5|) . while the only molecule for which the 
inclusion of anions has a significantly detrimental effect on agreement with observation is CsH. Both 
networks have difficulty reproducing the observed abundances of the sulfur-containing molecules, 
H2S, SO, C2S, SO2, and HCS + , due in part to the choice of t he low initial eleme ntal abundance for 
S + . A slightly higher choice can improve agreement strongly (jSmith et al.ll2004l ). With our current 
elemental abundance, SO, C2S, and SO2, all reach good agreement with observation at later times 
between 3 x 10 5 and 6 x 10 5 yr. Later times than the time of best agreement may pertain to clo udlets 
in TMC-1 other than the CP region, where most complex species peak (|Park et al.ll2006l ). The 
calculated abundances of several organic molecules are much lower than that observed, especially 
for CH3OH, CH3CCH, and CH2CHCN. All these species are difficult to synthesize using gas-phase 
processes only. It has been thought for some time that saturated species, in pa rticular, CH3OH, 
form on the surfaces of dust grains through the repeated hydrogenation of CO (jTielens &: Hagen 
1982 ). a claim backed up by some experimental evidence (jWatanabe et al.ll2004l ). 



The relatively small differences observed when anions are included in our models are partially 
due to their low abundance com pared with electron s , as shown in Figure [5] for Rate06. This 



is in contrast to the findings of IWakelam &: Herbstl (|2008l ) who compute that after ~ 10 2 yr, 



negatively charged small PAHs are the dominant negative charge carrier with an abundance almost 
two orders of magnitude larger than the electron abundance at steady state. This effect diminishes 
with increasing PAH size (and thus decreasing fractional abundance). With PAHs the dominant 
negative charge carriers, the overall ionization decreases because atomic positive ions are more 
readily neutralized by PAH~ than by bare electrons. This decrease results in increases in the 
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abundances of positively charged molecular ions, which in turn lead to increases in the abundances 
of unsaturated neutral species such as the cyanopolyynes and saturated species such as methanol. 
Despite the increase in the abundance of methanol, however, the computed abundance lies well 
below the observed one. Future models should investigate the effects of adding both negatively 
charged PAHs as well as molecular anions. 



VARIATIONS IN INITIAL ABUNDANCES 



In this work we have u sed oxygen-rich initial a b undances with a carbon - to-oxygen ratio of 



2002 



Smith et al 



2004 



~ 0.4. Previous mo dellers (jTerzieva fc Herbstl Il998l ; iRoberts fc Herbst 
Wakelam et al.ll20 06) have found that to get better agreement with observations of chemical abun 
dances in TMC-1 (C P), a higher rat i o is re quired, usually > 1. The most successful of the models 



listed is the work of Wakelam et al. ( 20061 ) who report that with a carbon-to-oxygen ratio of 1.2, 
86% of the abundances can be fit if the observations are given an uncertainty of a factor of five. 
These authors use as a criterion for success an overlap between Gaussian distributions of modeling 
and observational results for each species, the former obtained by a Monte Carlo analysis from 
uncertainties in both reaction rate coefficients and physical conditions. 

To test the effects that the addition of anions has on carbon-rich models of the chemistry of 
dark clouds, we have run models using both networks with a carbon-to-oxygen ratio of 1.2, with and 
without anions. Using Rate06 without anions, at the time of best agreement, 40 out of 51 species 
agree with obse r vation to within one order of magnitude, an improvement similar to that found by 
Wakelam et al. (j20od ). In contrast to the oxygen-rich initial conditions, we find that the addition 
of anions in the carbon-rich case worsens our agreement with observation for both networks. For 
Rate06 with anions, the agreement is reduced to 37 out of 53 species at the optimum time. The 
results with the OSU network are similar. In the carbon-rich case, the increased synthetic power 
obtained upon adding anions is magnified due to the increased carbon budget so that the calculated 
abundances of long-chain carbon molecules substantially overshoot the observed abundances, thus 
worsening the overall agreement with observation. Even in the case without anions, the better 
agreement obtained with carbon-rich abundances may be at least partially due to an edge effect, 
in which much of the carbon budget ends up in the longer-chained species. Since many of these 
species have not been observed in TMC-1, it is difficult to say whether the enhanced agreement is 
real or a result of the chemistry being limited to molecules with 11 carbon atoms or fewer. 



5. SUMMARY 

We have modelled the gas-phase chemistry of a dark cloud to investigate the role of molecular 
anions in the overall chemistry. We have used the simple pseudo-time-dependent gas-phase model, 
which can be criticized both for its lack of dynamics during the chemical evolution and for the 
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lack of processes involving dust grains. Nevertheless, our use of the simplified gas-phase model 
has enabled us to highlight the direct chemical effects of the addition of molecular anions on our 
existing standard networks. 

The inclusion of anions has only a slight positive effect on the success of both the Rate06 
and OSU networks in reproducing the observed abundances in TMC-1 (CP) with the pseudo-time- 
dependent model. The most significant effect of the addition of anionic species is the enhancement 
in the abundances of several families of carbon-chain molecules, namely, C n , C n H, C n H2, C n N, 
and HC n N, by providing new formation routes to the larger members of each family via smaller 
molecules. We get much better agreement with abundances observed in TMC-1 (CP) for the 
cyanopolyyne molecules, in particular, and would urge further experimental measurements of the 
rate coefficient and product channels for reaction ([5]), which we have shown to be an important 
formation route to the HC n N molecules. This 'catalytic' effect could be at least partially artificial 
because our chemistry is capped at molecules with 11 or fewer carbon atoms, so that we may 
be witnessing an 'edge' effect whereby molecules move up the chain through repeated addition 
of a carbon atom with the larger molecules becoming a 'sink' for carbon. To check if this is the 
cas e, it may be necess ary to add larger carbon-chain species to the chemistry, similar to that done 



by iMillar et al.l (|2000l ) for their chemical model of the shell of IRC+10216. Credence is given to 
our results, however, as the peak abundances of the molecules remain inversely proportional to 
molecular size even given this catalytic effect. 

We have also investigated the effect of using carbon-rich elemental abundances. Running 
models using both networks, with and without anions, with a carbon-to-oxygen ratio of 1.2, we 
find that we get improved agreement for TMC-1 (CP) over the oxygen-rich case for the carbon-rich 
models without anions. The inclusion of anions, however, worsens our agreement with observation 
due to the excessively enhanced synthesis of long-chain carbon species. 

Many of the large number of reactions involving anions that have been added to the reaction 
sets remain poorly understood. Although our results here show that the inclusion of anions does 
not have a dramatic effect on the chemistry, these conclusions must remain tentative pending 
further experimental work to determine more accurate rate coefficients. It is anticipated that 
in the near future, further anionic species will be discovered as laboratory measurements of the 
rotational spectra of candidate molecules increase in number, and we would expect that existing 
chemical networks will be expanded to incorporate molecular anions as standard. The current 
(anionic and non-anionic) releases of both Rate06 and the OSU network are available online at 
http: / / www.udfa.net /| and |http : / / www .physics . ohio-st ate, edu / ~eric/| respectively. 
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Fig. 1. — The abundances with respect to H2 density of (a) C n H _ for n = 4, 6, 8, and 10, (b) 
C~ for n = 4, 6, 8, 10, and (c) C n N~ for n = 3, 5 and their corresponding neutral analogs, as a 
function of time. 
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Fig. 2. — The anion-to-neutral ratio C n H~/C n H for n = 4, 6, 8, and 10, as a function of time. The 
dotted-dashed horizontal lines are the observed ratios for CqH~ and C§H~. 
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Fig. 3. — The abundances with respect to H2 density of (a) C n for n = 4, 6, 8, and 10, (b) C n H for 
n = 4, 6, 8, and 10, and (c) C n N for n = 3, 5, 7, and 9, as a function of time, with (dotted-dashed 
lines) and without (solid lines) anions. 
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Fig. 4. — The abundances with respect to H2 density of (a) C n H.2 for n = 4, 6, 8, and 10, and (b) 
HC n N for n = 3, 5, 7, and 9, as a function of time, with (dotted-dashed lines) and without (solid 
lines) anions. The dotted lines in the HC n N plot are results from the model using the larger rate 
coefficients for C n H~ + N -> HC n N + e". 
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Fig. 5. — The electron and total anionic abundances with respect to H2 as functions of time 
computed with Rate06. 
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Table 1. Initial Species and Elemental Abundances With Respect to Hydrogen Nuclear Density 



Element Initial Fractional Abundance 



H 2 


5.00 x 10 


He 


1.40 x 10 


C 


7.30 x 10 


N 


2.14 x 10 





1.76 x 10 


s+ 


2.00 x 10 


Si+ 


3.00 x 10 


Na+ 


3.00 x 10 


Mg+ 


3.00 x 10 


Fe+ 


3.00 x 10 


P+ 


3.00 x 10 


F 


2.00 x 10 


CI 


3.00 x 10 



References. - 



Graedel et al 



41982) 
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Table 2. Calculated Abundances of Families of Species With Respect to H 2 and Comparison 

with Observed Species in TMC-1 (CP) 



Rate06 OSU 
Species Observed 8, Without Anions With Anions Without Anions With Anions 



Anions 



C 4 H- 


< 2.3(-12) b 




l.l(-lO) 




l.l(-lO) 


C 6 H" 


1.2(-ll) c 




l.O(-lO) 




l.l(-lO) 


C 8 H- 


2.1(-12) c 




2.5(-ll) 




2.5(-ll) 


C 3 N^ 


< 7.0(-ll) b 




1.6(-12) 




3.6(-12) 


C„ 


c 4 




6.1(-09) 


3.8(-09) 


1.8(-08) 


1.6(-08) 


c 6 




6.7(-ll) 


3.2(-10) 


l.l(-lO) 


4.7(-10) 


c 8 




1.4(-11) 


1.4(-10) 


2.7(-ll) 


2.0(-10) 


Cio 




8.6(-13) 


7.9(-ll) 


4.3(-12) 


9.2(-ll) 


C„H 


C 4 H 


6.1(-08) c 


1.2(-08) 


1.2(-08) 


6.5(-09) 


9.8(-09) 


C 6 H 


7.5(-10) c 


2.1(-10) 


1.9(-09) 


1.8(-10) 


2.1(-09) 


C 8 H 


4.6(-ll) c 


1.3(-11) 


6.3(-10) 


1.3(-11) 


6.3(-10) 


CiqH 




1.7(-13) 


1.9(-10) 


1.5(-12) 


1.8(-10) 


C n H2 


C4H2 




3.2(-10) 


1.2(-09) 


1.4(-10) 


1.1 (-09) 


C6H2 




3.6(-ll) 


4.6(-10) 


3.6(-ll) 


5.9(-10) 


C 8 H2 




6.1(-13) 


1.6(-10) 


6.9(-13) 


2.0(-10) 


C10H2 




4.3(-17) 


3.2(-ll) 


3.6(-17) 


3.5(-ll) 
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Table 2 — Continued 



Rate06 OSU 
Species Observed 11 Without Anions With Anions Without Anions With Anions 



C n N 



C 3 N 


6.0(-10) b 


8.5(-10) 


l.l(-09) 


6.0(-10) 


1.0(-09) 


C 5 N 




3.8(-ll) 


9.7(-ll) 


3.3(-ll) 


1.5(-10) 


C 7 N 




4.4(-12) 


5.0(-ll) 


4.4(-12) 


6.2(-ll) 


CgN 




1.8(-13) 


1.7(-11) 


2.9(-13) 


2.0(-ll) 


HC n N 


HNC 3 


6.0(-ll) 


2.8(-ll) 


2.4(-ll) 


1.2(-10) 


1.5(-10) 


HC 5 N 


4.0(-09) 


3.6(-10) 


4.6(-09) 


2.2(-10) 


3.1 (-09) 


HC 7 N 


l.l(-09) 


5.4(-ll) 


2.1(-10) 


1.6(-11) 


1.2(-10) 


HCgN 


4.5(-10) 


3.9(-12) 


6.8(-ll) 


1.9(-12) 


3.6(-ll) 






Remaining 


Species Observed in 


TMC-1 




CH 


2.0(-08) 


5.5(-09) 


6.2(-09) 


6.6(-09) 


7.4(-09) 


NH 3 


2.0(-08) 


1.5(-08) 


1.3(-08) 


1.4(-08) 


1.3(-08) 


OH 


2.0(-07) 


4.0(-08) 


3.8(-08) 


2.7(-08) 


2.7(-08) 


H 2 


s^7.0(-08) 


1.2(-06) 


1.4 (-06) 


1.2(-06) 


1.4(-06) 


c 2 


5.0(-08) 


2.3(-09) 


2.9(-09) 


3.7(-09) 


3.8(-09) 


C 2 H 


2.0(-08) 


1.4(-08) 


1.6(-08) 


l.l(-08) 


1.6(-08) 


CN 


5.0(-09) 


2.6(-08) 


2.8(-08) 


8.3(-09) 


9.9(-09) 


HNC 


2.0(-08) 


2.4(-08) 


3.3(-08) 


2.0(-08) 


2.6(-08) 


HCN 


2.0(-08) 


3.4(-08) 


4.8(-08) 


1.9(-08) 


2.5(-08) 


HCNH+ 


2.0(-09) 


1.9(-10) 


4.5(-10) 


2.4(-10) 


4.7(-10) 


CO 


8.0(-05) 


1.0(-04) 


9.4(-05) 


9.6(-05) 


9.8(-05) 


N 2 H+ 


4.0(-10) 


l.l(-lO) 


l.l(-lO) 


1.9(-10) 


1.7(-10) 


HCO+ 


8.0(-09) 


4.0(-09) 


4.7(-09) 


4.0(-09) 


5.0(-09) 


NO 


3.0(-08) 


3.4(-08) 


2.4(-08) 


9.1(-09) 


8.9(-09) 


H 2 CO 


5.0(-08) 


l.l(-09) 


1.3(-09) 


8.6(-08) 


8.2(-08) 
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Table 2 — Continued 



Rate06 OSU 



Species 


Observed 3, 


Without Anions 


With Anions 


Without Anions 


With Anions 


CH3OH 


3.0( 


-09) 


5.1(-13) 


4.9(-13) 


1.2(-13) 


1.6(-13) 


2 


<7.7( 


-08) 


1.1 (-07) 


7.5(-08) 


7.3(-08) 


6.9(-08) 


H 2 S 


5.0( 


-10) 


7.2(-12) 


6.9(-12) 


9.2(-12) 


9.1(-12) 


C 3 H 


1.0( 


-08) 


1.4(-08) 


1.7(-08) 


1.0(-08) 


1.7(-08) 


C3H2 


1.1( 


-08) 


2.3(-08) 


1.5(-08) 


1.2(-08) 


1.1 (-08) 


CH 2 CN 


5.0( 


-09) 


2.4(-10) 


4.0(-10) 


2.2(-10) 


3.3(-10) 


c 2 o 


6.0( 


-11) 


4.3(-ll) 


7.5(-ll) 


1.5(-11) 


1.9(-11) 


CH3CCH 


6.0( 


-09) 


4.1(-12) 


4.2(-12) 


4.0(-12) 


4.6(-12) 


CH3CN 


6.0( 


-10) 


9.3(-10) 


1.9(-09) 


4.9(-10) 


7.9(-10) 


H 2 CCO 


6.0( 


-10) 


7.1(-09) 


9.4(-09) 


2.1(-08) 


2.3(-08) 


cs 


4.0( 


-09) 


4.9(-09) 


3.9(-09) 


7.1(-10) 


6.2(-10) 


/"ITT H T T /" \ 

CH3CHO 


6.0( 


-10) 


4.2(-ll) 


4.4(-ll) 


1.1(-12) 


1.3(-12) 


HCS+ 


4.0( 


-10) 


2.1(-12) 


2.9(-12) 


8.2(-13) 


9.8(-13) 


H 2 CS 


7.0( 


-10) 


1.8(-09) 


1.6(-09) 


2.4(-10) 


2.1(-10) 


r t /-' ~ \ /' \ /■ \ t T 

HCOOH 


2.0( 


-10) 


2.0(-09) 


2.4(-09) 


7.5(-10) 


1.0(-09) 


SO 


on/ 


-uy ) 


5.4(-ll) 


3.8(-ll) 


6.6(-ll) 


5.7(-ll) 


HC 2 NC 


5.0( 


-10) 






2.2(-10) 


3.8(-10) 


HNC 3 


6.0( 


-11) 


2.8(-ll) 


2.4(-ll) 


1.2(-10) 


1.5(-10) 


HC 3 NH+ 


1.0( 


-10) 


l.l(-ll) 


2.3(-ll) 


5.1(-11) 


1.3(-10) 


C 3 


1.0( 


-10) 


6.0(-10) 


7.5(-10) 


1.5 (-09) 


2.1 (-09) 


CH 2 CHCN 


4.0( 


-09) 


5.8(-13) 


1.3(-12) 


1.2(-13) 


2.0(-13) 


c 2 s 


8.0( 


-09) 


7.9(-10) 


7.6(-10) 


4.2(-10) 


5.0(-10) 


ocs 


2.0( 


-09) 


4.8(-10) 


4.6(-10) 


l.l(-lO) 


1.2(-10) 


so 2 


1.0( 


-09) 


5.4(-13) 


2.8(-13) 


8.2(-13) 


6.7(-13) 


CH3C4H 


4.0( 


-10) 


2.3(-ll) 


3.4(-ll) 


3.8(-ll) 


6.3(-ll) 


CH3C3N 


8.0( 


-11) 


3.0(-10) 


4.7(-10) 


9.8(-12) 


2.2(-ll) 


c 3 s 


1.0( 


-09) 


4.2(-10) 


4.0(-10) 


1.3(-10) 


1.6(-10) 


Agreement 






32/51 


35/53 


33/52 


36/54 



a As collated by ISmith et al.1 (|2004l ) unless otherwise stated 
^Thaddeus et al.1 feoo4 ) 
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Briinken et all (|2007|) 



Note. — a(b) represents a x 10 b 



Note. — Calculated abundances refer to the time of best agreement for each model, which is 
1-2 xlO 5 yr 



Note. — Normal type refers to modelled abundances which agree with observation to within 
one order of magnitude, bold type refers to those more than an order of magnitude smaller, and 
bold italic type refers to those which are more than one order of magnitude larger. 



